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We evaluate the transverse momentum spectrum ol J/ip (up to semi-hard momenta) in 
pA and AA collisions taking into account only the initial state effects, but resumming them 
to all orders in a^A 1 / 3 . In our previous papers we noticed that cold nuclear matter effects 
alone could not explain the experimental data on rapidity and centrality dependencies of the 
J/ip yield in AA collisions indicating the existence of an additional suppression mechanism. 
Our present calculations indicate that the discrepancy persists and even increases at semi- 
hard transverse momenta, implying a significant final state effect on J/ip production in this 
kincmatical domain. The QCD dipolc model we employ is only marginally applicable for J/ip 
production at mid-rapidity at RHIC energies but its use is justified in the forward rapidity 
region. At LHC energies we can quantitatively evaluate the magnitude of cold nuclear 
matter effects in the entire kincmatical region of interest. We present our calculations of 
J/ip transverse momentum spectra in pA and AA collisions at LHC and RHIC energies. 



I. INTRODUCTION 



Understanding the nuclear modification of J/ip production is a long-standing problem of high 
energy nuclear physics. The J/ip suppression was proposed as a litmus test of the quark-gluon 
plasma (QGP) formation [1] in the early days of relativistic heavy ion program. Since then it was 
realized that even though J/ip does offer a very valuable insight into the properties of the medium, 
the interpretation of the "anomalous" J/ip suppression is not straightforward. It is still not clear 
whether J/ip "melts" at the critical temperature [2, 3]. Moreover, J/ip production is affected by 
the "cold nuclear matter" (CNM) effects that tend to suppress it in a way that may resemble 
the suppression in QGP. The nature of these CNM effects is a subject of a controversy; however 
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it is clear that the coherence in the longitudinal direction must play a pivotal role. Since the 
knowledge of the CNM effects is crucial for quantifying the impact of the QGP on J/ip production, 
we have been motivated to perform the investigation of the CNM effects in a series of papers [4-7]. 
There we calculated the total cross section for production of J/ip in pA and A A collisions in the 
framework of the QCD dipole model [8]. Let us briefly explain the essence of this approach to J/ip 
production. At high energies, the cross section is dominated by the t-channel gluon exchanges. The 
properties of the J/ip wave function with respect to C and P transformations require that only an 
odd number of gluons are attached to the c and c pair making up the J/ip. Taking into account 
the large coherence length, one can realize that the dominant contribution to the J/ip cross section 
comes from multiple gluon exchanges along its path through the nucleus. In a large nucleus such 
that a 2 s A 1 /^ ~ 1 the gluon exchanges with different nucleons are parametrically enhanced [9-11]. 

The study of the total cross section as a function of rapidity and centrality in [7] indicated that 
the mechanism of the nuclear modification in CNM that we proposed could explain the experimental 
data on J/ip production in dA collisions at the RHIC. In A A collisions we found that it was 
responsible for a significant suppression of J/ip production, although it could not account for the 
entire effect. The total cross section is dominated by J/Y>'s produced at low transverse momentum 
due to the rapid fall-off of the differential cross section. A more detailed information is available 
through the study of the J/ip transverse momentum spectra that we perform in the present work. 
To understand the general form of the J/ip spectrum at semi- hard momenta and at high energies 
note that if the saturation momentum Q s is much larger than the quark mass, the geometric 
scaling [20-25] ensures that the cross section at a given centrality scales as da/d 2 p±d 2 b± ~ Qi/p±- 
The factor (Q 2 ) 2 indicates that at least two nucleons take part in the scattering. At high p± the 
coherence is lost, the geometric scaling breaks down and the power law of the p± dependence is 
determined by the relevant non-perturbative matrix elements instead of Q 2 . The transition from 
the semi-hard to the hard regime occurs at transverse momenta of the order of Qg/m 2 [22-25], but 
the detailed theoretical and phenomenological information about this region is scarce. The dipole 
model is applicable in the semi-hard region that we investigate in this paper. 

The paper is structured as follows. In Sec. II we derive the formulae for the J/ip transverse 
momentum spectrum in the semi-hard region by summing up all multiple interactions of the cc 
dipole with the nucleus compatible with the J/ip quantum numbers; we follow our analysis in 
[7]. Since we are not interested in the spectrum per se, but rather in its modification in the 
CNM, we re- write the cross sections for pA and AA collisions as a convolution of the cross section 
in pp collisions and the corresponding nuclear-dependent scattering factors. The main result is 
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given by (36), (37), (35), (31). In Sec. Ill we discuss the quasi-classical approximation and generalize 
our results to include the low-x evolution. In particular we found that the peculiar dependence 
of the scattering amplitude on the longitudinal coordinate makes it possible to derive a simple 
approximate formula for inclusive cross section given by (49), (52). In Sec. IV we perform the 
numerical calculations of the nuclear modification factor (NMF) using the DHJ model [12] for the 
forward dipole-nucleus scattering amplitude. Our results are displayed in Fig. 4, Fig. 5 and Fig. 6. 
Our calculation is in a reasonable agreement with the experimental data on dAu collisions at the 
RHIC. However, it tends to overestimate the NMF in AuAu collisions, especially at high pr- This 
probably indicates that the final state effects responsible for J/tp suppression grow at higher p?- 
One possible scenario for an additional nuclear suppression is due to the absorption in the quark- 
gluon plasma. An additional contribution to the nuclear suppression was recently investigated by 
one of us in [26, 27] where it was found that a strong magnetic field [28, 29] created in heavy ion 
collisions strongly suppresses high px J/ip's. It must be stressed that the dipole model can only 
be used for a qualitative estimate at mid-rapidity at RHIC energies because of the short coherence 
length; this is illustrated in Fig. 8 where we show the nuclear longitudinal formfactor. Neverthless, 
since the coherence length grows exponentially with rapidity, already at rapidity y = 1.7 and 
Pt < 2 GeV at RHIC the dipole model provides a fair estimate of the NMF. At the LHC the dipole 
model is applicable even at higher pt's. Therefore, a definitive test of the predicted enhancement 
of the final state effects at higher px can be done at the LHC. To this end, we present in Fig. 6 our 
calculation for the center-of-mass energy 7 TeV for pA and PbPb collisions. This is our prediction 
for the CNM effect at the LHC. 

II. J/ip PRODUCTION CROSS SECTION 

We separate the contributions to the cross section into three classes [7]: the scatterings that 
occur before and after the last inelastic interaction and the inelastic interaction itself. 

A. Multiple scattering before and after the last inelastic interaction 

The sum of diagrams depicted on Fig. 1 gives the contribution of an inelastic scattering of color 
octet before the last inelastic interaction. We have 

\ (V( Xl - x[) - V{0)) + \ (V(x 2 - x' 2 ) - V(0)) . (1) 





FIG. 1: One of the interactions before the last inelastic scattering. The diagrams that are complex conjugate 
to the first row of diagrams are not shown. The vertical dashed line denotes the last inelastic interaction 
when the cc pair is converted into the color-singlet state. The vertical solid line is the cut corresponding to 
the final state. 
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Multiplying (1) by pT(B)a, where p is the nuclear density, a is the total dipole-nucleon cross 
section, T(B) is the nuclear thickness, we get 
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where Q 2 is the gluon saturation scale given by (21) and p an infrared cutoff. Eq. (5) is the lowest 
order expansion of the dipole-nucleus scattering amplitude. It provides the initial condition for the 
low- a; evolution [10, 19]. This evolution erases the dependence of the scattering amplitude on the 
infrared scale p so that the amplitude becomes dependent only on the saturation scale - the effect 
known as the geometric scaling [20-24]. This allows to drop the logarithmic factors, as suggested 
by Golec-Biernat and Wusthoff [16] , and to write 



- l -Q 2 s (B) l( Xl - x ' 1 f + ( X2 - x ' 2 
Our notation is x 2 = x 2 = x 2 , . 
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FIG. 2: One of the interactions after the last inelastic interactions. Complex conjugate diagrams are not 
shown. 

Sum of the diagrams on Fig. 2 represents the contribution of an elastic scattering of color singlet 
after the last inelastic scattering. The corresponding factor is given by 

-Iq 2 s (B) [(»! - x 2 f + (x[ - x' 2 f] = -\q 2 s {B){t 2 + r' 2 ) , (7) 

where we again dropped the logarithmic factors. 

B. The last inelastic interaction 




FIG. 3: The last inelastic interaction. 

Sum of the diagrams on Fig. 3 represents the contribution of the last inelastic scattering con- 
verting the color octet into the color singlet. Neglecting logarithms we get 

- l -Q 2 s {B) [(x 2 - x[) 2 + (a* - x' 2 f - ( Xl - x[) 2 - (x 2 - x' 2 f] = -\q 2 s (B) rr>. (8) 
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This formula correctly describes the behavior of the total cross section, which receives its main 
contribution from the small-px region. However, since we are interested in the pj_-spectrum, (8) 
must be corrected by the logarithmic factors. Let us study the diagrams of Fig. 3 in more detail. 
The amplitude for J/ip production can be written as 

A(p±) = \?{pi_) , (9) 
P± 

where J 7 (p±) is the form- factor given by 

HP±) = J *G(k) V v (k -p) = 2vra s £ dz j ^-<S>(r, z) (e'^ - e 1 ^) . (10) 

Here an d are the wave functions of the s-channel gluon and J/ip correspondingly, and $ 
is explicitly given by (28). In deriving (10) we assumed that the coherence length is much larger 
than the nuclear radius. This is strictly true in the Regge limit. However, at RHIC and LHC the 
effects of a finite coherence length may be important, especially for heavy particles and/or high 
transverse momenta. In Appendix A we compute the longitudinal form factor Tl that is a measure 
of the longitudinal coherence, and discuss the kinematics in which our results in this paper are 
applicable. 

Introduce the impact parameter representation of the scattering amplitude 

A(b) = I W? A ^ l= Yj T- 1 ^*(r)(e-^'-e->>>-ir><) . (11) 
Then the cross section for scattering on a given nucleon reads 

= \A( Pi _)\ 2 = j d 2 bj d 2 b' e l ^ b - b '^A{b)A*{b') , (12) 

To obtain the cross section for scattering on a nucleus we need to average the forward scattering 
amplitude over all possible positions of the nucleon in the nucleus. We have 



d 2 p 

where 



d ^ = j ' d 2 bj d 2 b' e^O'V (A(b)A*(b')) . (13) 



(...) = J d 2 b a P T(b a )A(b - b a )A*(b' - b a ) . (14) 

Here b a is the transverse position of the nucleon with respect to the center of the nucleus, p is the 
nuclear density and T(b) is the nuclear thickness. For a big nucleus pT(b a ) « A/Sa is a constant, 
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with Sa being the nuclear cross-sectional area. We have 

„2 r r r Alt, 



(A(b)A*(b>)) =^7^2 j dT j ^ (e-^HrH _ e -^-6 -|r.).^ 



^ / i(V-6„+ir')-/' _ e i(6'-& 
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We introduced a convenient notation 



i dT= j! dz i s* (r ' z) r ^ / ? ^ z ° (i7) 

Substitution of (16) into (13) gives 



If we were interested in the total cross section, then integrating (18) over p± would give 
do g A 



d 2 b 



aiA-JdrJ^ (e^^'H + e-i^-r'yi _ e -ifrW)-l _ ^(t^'K) (19) 

o A r m r d 2 £ l . 



where we expanded the exponents in (19) because only small momenta I ~ 1/Ra give the loga- 
rithmic contribution to the integral. Note that formally, (20) averages to zero. This is because 
J ftp is C and P odd and therefore there must be an odd number of gluons hooked to the fermion 
line, whereas (20) accounts only for two gluons. However, once the multiple scatterings are taken 
into account, the r ■ r' term selects the terms with the required parity form the scattering factors 
before and after the last inelastic interaction as we demonstrate below. Thus, we think of (20) as 
the projection operator onto the l - ~ state. Defining the saturation momentum 

Q 2 = ^a 2 sP T{B) (21) 

and dropping the \n(l/ fi\r ± r'\) factors we cast (20) into the form 

dcr g A 



j dTQ 2 s \r.r' . (22) 



d 2 b 

This is precisely the approximation in which our formula (8) is valid. We observe that in general, 
expansion (20) cannot be performed for the differential cross section (18). Indeed, in (18) at high 
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p±, A = b — b' is small and therefore I is large. Hence, we write (18) as 



dogA A a 



^2 j d 2 Ae ip -* J dTn J(r,r', A). 



(23) 



d 2 B d 2 P± S A (2 
where we introduced an auxiliary function J as 

J(r, r', A) = J ^-e'^ e (e^ r ~ r '^ + e^'^ - e'^^ 1 - jW+V)-^ (24) 
Using formula (B4) of Appendix B we find 



J(r,r',A) ( 



I 2 

2 



) + A| 

+ Q (r + r - ) + A )^ln ^L^ 

+ (l(, + r 0-A) 2 m ' 



(25) 



I 2 

In particular, J(r,r', 0) = r • r' modulo the logarithmic factors. 

C. Summing up the multiple scatterings 

Suppose that the longitudinal coordinate of the last inelastic scattering is Then using 
(6), (7), (23) and summing over all possible number of interactions, we find that the scattering 
amplitude for J/tp production in pA collisions is given by 

T pA ^„x (r, r>, B, A) = J™ \ f J(r , r> , A) e & (-T]| e -§Q^)(i-|) 

4J(r,r',A) f i Q2(B)[(r-r') 2 +4A2] _ -§ Oj(B)(r 2 +r' a ) \ ( 2 6) 
( r + r /)2_ 4A 2 \ e e J' ^ 

The cross section for inclusive J/Y> production in pA collisions reads 



d&pA^J/ipX 



where 



x | ^c^r^j/^^^A,^) (27) 



$(r,z) = ^-|= {m^ o (m c r)0 T (r,z) - [z 2 + (1 - z) 2 ] m c ^i(m c r)^.0r(r,^)} (28) 
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with 



f>T(r, z) = Ntz(1 — z) exp 



(29) 



2Rj^ _ 

and where Nt = 1-23, BJf, = 6.5 GeV -2 . Integrating over p± we reproduce the result of [7]. 

It is well-known that the leading order calculation does not properly describe the transverse 
momentum distribution of hadrons. Moreover, the high-px tail of such distribution stems from 
the short distance "hard" processes, which violate the geometric scaling. These hard processes 
factorize from the rest of the scattering amplitude in a number of recently studied reactions [30]. 
Since we are primarily interested in the nuclear modification of J/ip production in pA and AA 
collisions as compared to pp ones, we propose a model for the function J of (24), which retains the 
essential features of (25) and additionally describes the J/tp spectrum in pp collisions. In view of 
(20) we suggest to write 

J(r,r',A)wr-r'F(A), (30) 



F(A) = 1 / W(Pi) e -,, V|)1 j (31) 

<T vv^MX J « P± 



where the short distance effects are described by the function F(A), which is defined as 

with Upp being the total inelastic pp cross section. The factor r ■ r' in (30) selects (after integration 
over the angle between r and r') only those contributions to the cross section that contain odd 
number of gluons hooked up to the fermion line in agreement with the J ftp quantum numbers. 
According to our proposal we model the scattering amplitude (26) by 

T pA ^Mry,B, A) - (Xyiti* { e -" e; ' B "<"- ,) ' +4A ' 1 - .-WW*")} . (32) 

Expanding (32) to the leading order in Q 2 S we obtain the scattering factor for pp — > J/ipX process: 

T pp ^, WX (r, r', B, A) =\r • r'Q 2 s F(A) + • r'(-3r 2 - 3r' 2 + 2r • r' - 4A 2 )Q^ (A ) . (33) 

Upon averaging over the angle between the dipoles r and r' we get 

(T pp ^ x {ry,B,A)) = l(r ■ r') 2 Q 4 s F(A) . (34) 

Appearance of Q$ indicates that at least two scatterings are required to produce J/ip in accordance 
with the C and P odd nature of the J/tp wave function. In our model F(A) completely determines 
the pj_-spectrum.* Our goal is to calculate the modification of the J/ip spectrum in the cold nuclear 
medium. 



Integration with the complete function J selects dipole sizes r, r' ~ l/p± so that the cross section falls off as Qi/p± 
at large p± . 
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Generalization of the scattering amplitude to AA collisions reads 
T AlA ^j m (r,r ,B, A) - ^ ^ + ^ (r + r / )2 _ 4A ^ 

X | e -^(Q?l+Q? 2 )[(»-'-') 2 +4A 2 ] _ e -I(Q2 1+ Q2 2 )( r 2+ r '2) j _ (35) 

The corresponding cross section for inclusive J/tp production in AA collisions is given by 

^WB 2 =[ iZ I *« r '*£ ^^Mry, B, P± ), (36) 

where 

T Al A^j/*x(ry,B,p ± ) = J ^2e iA -PT AlA2 ^ jm (r,r',B,A) (37) 
We will simplify the nuclear profiles by the step-function. 

III. QUASI-CLASSICAL APPROXIMATION AND EVOLUTION EFFECTS 

A. Quasi-classical approximation 

Numerical integration over A in (37) is difficult as the integrand is an oscillating function. In the 
quasi classical approximation we can avoid this problem if we neglect the logarithmic dependance of 
the saturation momentum on the dipole size as discussed in the previous section. Let us introduce 
the following Fourier representation 

4r -/ L-^QUBMr-rnW] _ e -|<K(B)(^) j = f ^ ^A-k R{ j g fc) 

(r + r'Y — AA Z I ) J 

(38) 

Using (35), (37) and (38) we cast the scattering amplitude in the form of a convolution 

T AlMm (r, r>, B,p)= ° F 2 ®} X f% 2 f F k H(r,r',B,p-k)d 2 k ± , (39) 

where 

Fk = J |^ elA '" F(A) • (40) 

Note, that by (31) / d 2 kF k = F(0) = 1. The A-integral in the inverse of (38) can be done analyt- 
ically. Let us invert (38) and use (26) to write it as an integral over the longitudinal coordinate 
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t' = £/T{B): 

d 2 A ^. kr r , (Q 2 sl + Q 2 s2 ) 



H < r >'> B >*=Jwr~ ■ ■ is 

X C ^' e -K^i+^2)(A 2 +i(r-r-') 2 )r e -|(Q 2 1 +Q 2 2 )(r 2 +r' 2 )(l- 5 ') ^ 

Jo 

Jo 16tt£' V ; 



Finally, the cross section reads 



d 2 p ± dyd 2 B 1 d?B 2 



XJ ( 2 2jUfc± 2 e -|(Q? 1 +Q 2 2 )|(--r-') 2 «' e -|(Q? 1 +Q 2 2 )(r 2 +r'2)(l- 5 ') (43) 

where Jo is a modified Bessel function and B\ 2 are the impact parameters of the two nuclei. 

B. Low- a; evolution 

Explicit integration over A is not possible if the amplitude is evolved at low x. Nevertheless, 
we found a very good approximation to our formulas even in this case. Let us again re-write the 
cross section (36), (35) as an integral over the longitudinal coordinate £' using (26). 

where we introduced an auxiliary function S as 

m=l 1 dzJ d ^r, Z )£d Z 'J d ^(r', Z ')r.r' 

It is not difficult to see that E varies with £' much slower than the exponential factors in (44). It 
starts at 2(0) = (note the integration over the angle between r and r') and rises at most as a 
power towards a constant 5(1). We therefore approximate (44) as 



do'A 1 A 2 ~^J/ipx ^ Cf 
d 2 p ± dyd 2 B 1 d 2 B 2 



F - 2 Q 2 sl Q 2 s2 (s> jf 1 d? J |^ e ^- PF(A)e -i(e^ 2 )A 2 r t 
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where the average value of S is given by 



Jo 



( 

The pj_-integrated cross section reads using (46) 

d (7 A 1 A 2 ^J/i>X _ Cf 



(47) 



dyd 2 B 1 d 2 B 2 4a s ir 2 
Now, employing (48) in (46) we obtain 



(48) 



da 



A 1 A 2 ^J/ipX 



da 



AiA 2 ^J/ipX 



-1 [^ e ^F(A)±(l- 

\ l + QlJ (2.)^ *^A 2 V 



-i(Q2 1+Q 2 2 ) A 2 



e 4 



d 2 p±dyd 2 B 1 d 2 B 2 dyd 2 B 1 d 2 B 2 Q 2 sl 
We can incorporate the low- a; evolution effects by the following replacements 



) • (49) 



e-lQy 2 l-N F (r,b,y) 



e i6 



Therefore, the cross section reads 

d&A 1 A 2 ^J/ipX _d< J A 1 A 2 ^J/ipX 



l-N A (r/2,b,y) 



(50) 
(51) 



d 2 P± dyd 2 B 1 d 2 B 2 dyd 2 bd 2 B J (2vr) 



d 2 A 



P F(A) 



x 1 



l-NW^Buy)] [l-N® (A, B 2 , -!/)]) 
(l- [l-AT^(A,B 1)y )] [l-ivf (A,S 2 ,-y)])^ 



(52) 



where LT stands for the "leading twist" which is the leading term in small A expansion. The result 
for the pj_-integrated cross section can be found in our previous paper [7]. 

The experimental data is expressed in terms of the nuclear modification factor (NMF). It is 
defined as 

da a 



Ra x A 2 — 



1 2 dyd 2 p± 



(53) 



where S stands for the overall area of two colliding nuclei. Since the mechanism of J/tp production in 
pp collisions remains elusive, we follow our approach in the previous publications and approximate 



da, 



pp^rJ/ipX 



= C 



da 



AA-^J/tpX 



(54) 



A=l 



dy dy 

with C = const. We fix the constant to provide the best description of the pp and dA data. We 
found in the previous paper [7] that is close to unity, so we set C = 1 in the present paper as well. 
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IV. NUMERICAL CALCULATIONS 

Experimental data for the differential cross section of J ftp production in pp collisions can be 
parameterized in the following form [31]: 

<Jppd^p± V PoJ 

where M and po are positive constants and a pp is the total inelastic cross section. From (40) it 
follows that 

F(A)= i ^K 5 ( Po A). (56) 

We calculate the spectrum in pp — > J/i/jX process by setting A\ = A2 = 1 in (43). Although the 
saturation momenta are small in this case, there is still a small influence of higher twist terms on 
the final spectrum which results in a shift of the average transverse momentum, which according 
to (55) equals = Po/4. Because of this we cannot directly use the experimentally determined 
value of po, but rather fit the constants po and N to the data [31, 34]'. 

We performed the numerical calculations using the DHJ model of the dipole scattering amplitude 
[12]. The DHJ model is an improvement of the KKT model [13, 15] that takes into account the 
change in the anomalous dimension of the gluon distribution function due to the presence of the 
saturation boundary [14] and takes into account some higher order effect. It successfully describes 
the single inclusive hadron production in dA collisions in the relevant kinematic region. In this 
model, the dipole scattering amplitude is parameterized as follows 

V A (r,0,y) = l-exp{-^ (r 2 Q 2 ) 7 } . (57) 

The gluon saturation scale is given by 

Q 2 S = A 2 A 1 / 3 e Xy = 0.13 GeV 2 e Xy N coll . (58) 

where the parameters A = 0.6 GeV and A = 0.3 are fixed to the low- re DIS data [16] and are 
consistent [17, 18] with the RHIC and LHC results on hadron multiplicities. The anomalous 
dimension reads 

ln(M 2 /Q 2 ) , N 

7 = 7 + 1_ 7 1 .T 5 59 



t Since we are concerned with low pr J/tp's we fit po and M to the ALICE collaboration data [34]. 
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FIG. 4: (Color online). Nuclear modification factor for J/ip's vs p± in GeV at y/s = 200 GcV, y = in 
AuAu for centralities (a) 0-20%, (b) 20-40%, (c) 40-60%, (d) 60-92% and in (e) minbias dAu. Data is from 
[32, 33]. 

where M± = ypj_ + 4m 2 , 7 S = 0.628 is implied by theoretical arguments [14] and d = 1.2 is 
fixed by fitting to the hadron production data in dA collisions at the RHIC. Y = ln(l/x), with 
x = me~ y / y/s. The quark dipole scattering amplitude is given by 



N F (r,0,y) = 1 - y/l-N A (r,0,y) 
In the DHJ model (52) reads 



(60) 



da 



A 1 A 2 ^J/tpX 



da 



AiAi->J/il>X 



d 2 A 



e iA - p F(A) 



d 2 p±dy dy J (2tt) 2 

The results of our numerical calculations using (61) are presented in Fig. 4, Fig. 5 and Fig. 6 
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FIG. 5: (Color online). Nuclear modification factor for J/ip's vs p± in GeV at y/s — 200 GeV, y = 1.7 in 
AuAu for centralities (a) 0-20%, (b) 20-40%, (c) 40-60%, (d) 60-92% and in (e) minbias dAu. Data is from 
[32, 33]. 



for the center-of-mass energies \fs = 200 GeV and \fs = 7 TeV. As mentioned in the Introduction 
they indicate that the final state effects on J/ip production increase with p±. 



V. SUMMARY 



In this paper we derived the p± spectrum of J/^'s produced in pA and AA collisions in the 
framework of the dipole model. We took into account the strong coherence effects in the cold 
nuclear medium, but entirely neglected the final state effects. We used a phenomenological model 
for the scattering amplitude to numerically investigate the transverse momentum dependence of 
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FIG. 6: (Color online). Nuclear modification factor for J/ip's vs p± in GeV at y/s = 7 TeV in PbPb for 
rapidities (a) y = and (b) y = 3.25. Each line corresponds to a different centrality bin; from bottom to 
top: 0-10% (solid red), 10-20% (dashed green), 20-30% (dash-dotted blue), 30-50% (solid purple), 50-80% 
(dashed magenta) and in minbias pPb (solid brown) . 



the nuclear modification factors. Our results provide a useful reference for evaluation of the contri- 
bution of the final state effects to the J/ip suppression; they are reasonable agreement with RHIC 
data on J/ip production in dAu collisions. The LHC pA data will be an important test of our 
approach based on the gluon saturation. We find that the J/ip suppression that originates from 
the initial state (cold nuclear matter) effects increases at the LHC energies compared to RHIC. 
Meanwhile, the experimental data on AA collisions indicate [35] that J/ip's are suppressed less at 
LHC than at RHIC. We have not found a solution for this problem - in fact, our results exacerbate 
it. 
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Appendix A: The longitudinal form factor Ti, 

The longitudinal form factor J^l is a measure of the coherence of a high energy process. It is a 
function of the longitudinal momentum transfer p z and is defined as 

?L(p z ) = \j d% l°° c%p(b,0e ip *t (Al) 
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where £ is the longitudinal coordinate. In the coherent regime p z £ <C 1 implying that J~l = 1. 
Denote the momenta of gluon, quark and antiquark as q = (cj,0,w), k\ = (zu>, k\, k\ z ) and &2 = 
((1 — z)u,p — fei, A;2 2 ), see Fig. Al. The longitudinal momentum transfer p z can be evaluated as 

fe? + m 2 (p — fci) 2 + m 2 , , . 

P* = & - fciz - fe, ~ + % n U , A2 

2zuj 2(1 — z)w 

The relative transverse momentum of c and c is small compared to the typical transverse momen- 




FIG. 7: A diagram contributing to the form factor. Dashed line represents all possible gluon exchanges. 

turn of the J/ip center of mass, viz. \k\ — k-z\ ~ a s m <C p±. Therefore, we can approximately set 
ki± ~ k2± ~ p±/2. Considering that typically z ~ 1/2 we derive 

Pz = ~ [p\ + (2m) 2 ] = ^e~y [p\ + (2m) 2 ] . (A3) 

Consider the "hard sphere" model of the nucleus: 

Pip, = ^ o(\t\- - 6 2 ) , (A4) 

where 6 is the step function. Then, substituting (A4),(A3) into (Al) we find that the longitudinal 
form factor reads 

3 

•F L = j- [sm(p z R A ) - p z R A cos(p z R A )} ■ (A5) 

The cross section is proportional to J-l which is plotted in Fig. 8. We see that the dipole model 
works well at the LHC and at rapidity y = 1.7, p± < 3 GeV at the RHIC, whereas at y = and 
the RHIC it gives a qualitative estimate at best. 

Appendix B: Evaluation of the integrals in (24) 

Consider a typical integral appearing in (24) 

= [ ^e«« (Bl) 
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FIG. 8: (Color online). Square of the longitudinal form factor as a function of transverse momentum of 
J/ip in GeV's. Lines from bottom to top: y/s = 0.2 TcV, y = (blue, solid), y/s = 0.2 TeV, y = 1.7 (red, 
dashed), y/s = 2.76 TeV, y = (green, dot-dashed), y/s — 5.5 TeV, y = (brown dashed). 

This integral is quadratically divergent in the infrared region. However, the quadratic divergent 
terms cancel between the four terms appearing in (24) as is evident in (20). Therefore, we are 
interested only in terms that diverge at most logarithmically. To find those, take Laplacian of I{x) 

dll{x) = -j ^ e -^ = - 27 rln-L (B2) 

Note that I{x) depends only on the absolute value of x. Using the polar coordinates we cast (B2) 
in the form 

if (B3) 

X OX \ OX J XfJ, 

Integrating this equation yields 

/(z) = -(2vr)^ln^- + ... (B4) 
where elapses indicate the divergent terms independent of x. 
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